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Abstract 

Atmospheric Cerenkov Technique is an established methodology to study TeV en- 
ergy gamma rays. However the challenging problem has always been the poor signal 
to noise ratio due to the presence of abundant cosmic rays. Several ingenious tech- 
niques have been employed to alleviate this problem, most of which are centred 
around the Cerenkov image characteristics. However there are not many techniques 
available for improving the signal to noise ratio of the data from wavefront sam- 
pling observations. One such possible technique is to use the Cerenkov photon arrival 
times and identify the species dependent characteristics in them. Here we carry out 
systematic monte carlo simulation studies of the timing information of Cerenkov 
photons at the observation level. We have parameterized the shape of the Cerenkov 
shower front as well as the pulse shapes in terms of experimentally measurable 
quantities. We demonstrate the sensitivity of the curvature of the shower front, 
pulse shape parameters as well as the photon arrival time jitter to primary species 
and show their efficiency in improving the signal to noise ratio. The effect of limiting 
the Cerenkov telescope opening angle by using a circular focal point mask, on the 
efficacy of the parameters has also been studied for each of the parameters. Radius 
of the shower front, pulse decay time and photon arrival time jitter have been found 
to be the most promising parameters which could be used to discriminate 7— ray 
events from the background. We also find that the efficiency of the first two param- 
eters increases with zenith angle and efficiency of pulse decay time decreases with 
increasing altitude of observation. 
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1 Introduction 



Atmospheric Cerenkov Technique (ACT), is a well established and a unique 
method for the astronomical investigation of Very High Energy (VHE, also 
referred to as TeV) 7— rays. It is based on the effective detection and study 
of the Cerenkov light emitted by the secondary particles produced in the 
extensive air showers initiated by the primary 7— ray [1-4]. 

As is well known, the abundant charged cosmic ray particles generate Cerenkov 
light akin to that produced by the 7— rays as a result of which the 7— ray 
signal is buried in a vast sea of cosmic ray background. This has been a major 
difficulty in applying the atmospheric Cerenkov technique successfully for 7— 
ray astronomy. 

The first generation Cerenkov telescopes addressed the problem of cosmic ray 
background by matching the telescope aperture to the angular size of the 
Cerenkov flash [5]. In addition, many of the potential astrophysical sources 
of VHE 7— rays are expected to produce modulated signals (e.g. 7— ray 
pulsars, X-ray binaries, cataclysmic variables etc). Thus techniques based on 
well established methods of phase sensitive detection were used to enhance 
the detection of 7— rays against the randomly arriving background [6,7]. 

This technique cannot be applied however for steady sources of 7— rays like the 
Blazars, normal galaxies or unidentified 7— ray sources discovered at GeV en- 
ergies. There are, however, differences in the detailed structure of the detected 
Cerenkov flash due to the fact that the cosmic ray background protons are 
spatially and temporally isotropic while the 7— rays emanate from the point 
sources. In addition, the physics of the hadronic cascades initiated by the pro- 
tons in the earth's atmosphere is different from the electromagnetic cascade 
generated by the 7— rays. Detailed monte carlo simulations have been em- 
ployed to study these differences. Techniques based on the shape of the lateral 
distribution of the Cerenkov light pool to delineate the source direction, have 
been studied [8,9] despite large fluctuations in the measured Cerenkov pho- 
ton densities [10,11]. There are several Atmospheric Cerenkov arrays designed 
precisely to apply these techniques to ground based VHE 7— ray astronomy 
[12-17]. On the other hand the imaging technique has been shown to be suc- 
cessful and also has been demonstrated to work reliably in detecting emission 
from several TeV 7— ray sources [18,19]. 

The spatial and temporal properties of the Cerenkov photons also contain 
valuable information on the development and propagation of the EAS in the 
atmosphere. As a result, systematic studies of these photons as received at the 
observation level could lead to the development of techniques to distinguish 
between hadronic or photon primaries. 7— ray primaries for example, develop 
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higher in the atmosphere resulting in compact images as well as narrower 
pulses. One would therefore expect correlations between different pulse shape 
parameters and imaging parameters. 

Extensive studies have already been carried out in this regard using detailed 
simulation techniques. Most of these studies were carried out at higher energies 
with the aim of studying the elemental composition of cosmic rays at these 
energies [20]. The possibility of applying pulse shape discrimination technique 
to improve the signal to noise ratio in the data from Crab nebula was first 
demonstrated by Turner et al. [21]. Their study was solely based on the pres- 
ence or absence of kinks and other anomalies like long trailing edges of the 
Cerenkov pulses indicating hadronic origin. These criteria were not based on 
any systematic simulation studies nor were they conveniently parameterized 
so that it could be used by others. 

Several characteristics of the hadron- & 7— ray showers which were suggested 
as possible discriminators in the past have been met with limited success. Some 
of these are the presence of penetrating particles among the EAS secondaries 
[22], the ultraviolet excess in hadron initiated showers due to their proximity 
to the observation level [23] and time duration [24,25]. On the other hand 
methods developed for improving the angular resolution technique [26] are 
still in use and are one of the important methods of discrimination for non- 
imaging arrays which are generally spread out. 

A brief description of the efficacy of the pulse shape discrimination in VHE 
7— ray astronomy was made by Patterson & Hillas [27]. Their study was 
based solely on the presence of structure on the leading or trailing edge of the 
Cerenkov pulses as suggested before [21]. However no systematic studies of 
the rise or decay times, FWHM etc were carried out. Recently Roberts et al. 
[28], developed a technique based on the temporal Cerenkov pulse shape and 
showed that the use of rise time at large zenith angles and FWHM at smaller 
zenith angles are effective discriminators. However these studies are carried 
out at larger zenith angle (> 35°) and hence the conclusions are relevant only 
at larger primary energies (> 40 TeV). On the other hand Cabot et al. [29] 
suggest possibility of identifying the muonic component, which precedes the 
electromagnetic component, of the Cerenkov pulse in order to separate the 
7— ray signal from the background. This could be a very efficient technique 
in experiments based on wavefront sampling technique. But once again this 
method works only at higher primary energies (> a few TeV) where a signif- 
icant number of muons are produced. 

In an earlier work[ll] functional fits have been carried out to the spherical 
shower front to demonstrate that the radii of curvature are equal to the height 
of the shower maximum at all observation levels. It has also been demonstrated 
in that paper that the Cerenkov photon arrival time distributions can be well 
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fitted by a lognormal function. Using these fits, radial and spectral variation 
of pulse shape parameters have been studied systematically. 

In the present work we plan to make a systematic study of the temporal 
and spatial profile of Cerenkov light from lower energy primaries both from 
pure electromagnetic cascades as well as hadronic cascades generated by TeV 
energy primaries. The question we are trying to answer is whether the observed 
differences in experimentally measurable temporal information could be used 
to separate electromagnetic component from the hadronic background. 

In §2 of this paper details of simulations are given, followed by definition of 
figure of merit for discrimination between 7— rays and cosmic rays in §3. In §4 
we discuss the results from the analysis of the curvature of the Cerenkov front 
whereas in §5 we present results based on the pulse shape parameters and 
identify the more sensitive of them. In §6 we present the detailed study of the 
photon arrival time jitter. The dependence of the parameters on the telescope 
opening angle, altitude of observation, incident angle of the primary as well 
as cosmic ray species and the shower core distance are discussed in §7. A brief 
discussion of the results is presented in §8 and conclusions are summarized in 
59. 



2 Simulations 

A package called CORSIKA (version 560), [30,31] has been used to simulate 
Cerenkov light emission in the earth's atmosphere by the secondaries of the 
extensive air showers generated by cosmic ray primaries or 7— rays. This pro- 
gram simulates interactions of nuclei, hadrons, muons, electrons and photons 
as well as decays of unstable secondaries in the atmosphere. It uses EGS4 code 
[32] for the electromagnetic component of the air shower simulation and dual 
parton model for the simulation of hadronic interactions at TeV energies. The 
Cerenkov radiation produced within the specified band width (300-650 nm) 
by the charged secondaries is propagated to the ground. The US standard 
atmosphere parameterized by Linsley [33] has been used. The position, an- 
gle, time (with respect to the first interaction) and production height of each 
photon hitting the detector on the observation level are recorded. 

In the present studies we have mainly used Pachmarhi (longitude: 78° 26' E, 
latitude: 22° 28' iV and altitude: 1075 m) as the observation level where an 
array of Cerenkov detectors each of areaQ 4.35 m 2 is deployed in the form of a 
rectangular array. We have assumed 17 detectors in the E-W direction with a 

1 This is the total reflective area of 7 parabolic mirrors of diameter 0.9 m deployed 
paraxially on a single equatorial mount. 
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separation of 25 m and 21 detectors in the N-S direction with a separation of 20 
m. This configuration, similar to the Pachmarhi Array of Cerenkov Telescopes 
(PACT) [12] but much larger, is chosen so that one can study the core dis- 
tance dependence of various observable parameters. Monoenergetic primaries 
consisting of 7— rays, protons and iron nuclei incident vertically on the top of 
the atmosphere with their cores at the centre of the array have been simulated 
in the present studies. The showers simulated in this study have a fixed core 
position which is chosen to be the detector at the centre of the array. The 
resulting Cerenkov pool is sampled by all the 357 detectors which are used to 
study the core distance dependence of the parameters studied here. All the 
telescopes are assumed to have their optic axes aligned vertically. 

An option of variable bunch size of the Cerenkov photons is available in the 
package which serves to reduce the requirement of hardware resources. How- 
ever since we are interested in the fluctuations of each of the estimated ob- 
servables, we have tracked single photons for each primary at all energies. 
Multiple scattering length for electrons and positrons is decided by the pa- 
rameter STEPFC in the EGS code which has been set to 0.1 in the present 
studies [34]. Wavelength dependent absorption of Cerenkov photons in the 
atmosphere is not however taken into account. The present conclusions are 
expected to be independent of photon wavelengths. 



3 Figure of merit of a parameter 

Figure of merit of a parameter that can distinguish between VHE 7— rays 
and cosmic ray hadrons depends primarily on two factors. Firstly, it should 
accept most of the 7— rays and secondly it should be able to reject most of the 
hadrons. In general, this figure of merit could be a function of primary energy. 
In the present work we define such a figure of merit which is often called as 
quality factor, as [28] : 




(1) 



where iVJ is the number of 7 rays accepted, 



Nj, is the total number of 7 rays, 



N£ r is the number of background cosmic rays accepted and 



Njf is the total number of background cosmic rays. 
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The quality factor thus defined is independent of the actual number of 7— rays 
and protons recorded. 

In this paper, whenever a pair of distributions of a parameter under study for a 
hadron and a 7— ray primary are shown the threshold value of the parameter 
is indicated by a vertical line. The threshold is chosen such that it yields 
maximum quality factor subject to the condition that the accepted fraction 
of 7— rays is > 30% and that of protons is > 1%. The errors on the qulaity 
factors shown in each case are statistical only. 



4 Shower front parameters 

It has been shown long ago that the radius of curvature of the Cerenkov light 
front is strongly correlated with the height of shower maximum from the ob- 
servation level [35]. This has been found to be true for different species of 
cosmic rays [11]. For photonic primaries the height of shower maximum is 
decided by the radiation length in the atmosphere while that for hadronic 
primaries it is decided by the interaction length which in turn depends on the 
interaction cross-section in air. Hence the radius of curvature could be species 
specific. Therefore we have investigated the possibility of using the fitted ra- 
dius of curvature of the spherical Cerenkov front as a parameter to distinguish 
between 7— ray and proton initiated showers. The details of a spherical fit to 
the mean arrival times of Cerenkov photons at detectors sampling the front at 
various core distances are described in detail in [11]. Fig. 1 shows the distri- 
bution of the fitted radii for different primary species of various energies. The 
primaries are incident vertically at the top of the atmosphere except wher- 
ever mentioned. Mono-energetic 7— ray and hadron primaries of comparable 
Cerenkov yield are simulated and compared. We also simulated showers, both 
for 7— rays and protons whose energies are selected randomly from a power 
law distribution of a differential slope of -2.65. This would also simulate the 
real data as one would record in an experiment. While the slope of the 7— ray 
spectrum could be different from what is chosen here, the quality factor does 
not depend on the number of showers and hence independent of the spectral 
slope except when it shows strong energy dependence. The energy bandwidth 
chosen here are 500 GeV - 10 TeV for 7- rays while it is 1 TeV - 20 TeV for 
protons. 

Quality factors, as defined above, have been estimated for the present sam- 
ple from the distributions of the estimated radii of curvature for 7— ray and 
hadronic primaries as shown in figure 1. These are given in table 1. The thresh- 
old values in each case are indicated as vertical lines in figure 1. Fraction of 
total number of 7— ray and proton showers that pass this cut, i.e., fraction 
of showers with radii below the threshold as well as the number of showers 
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Table 1 

Quality of radius of curvature of the spherical photon front as a discriminating 
parameter for vertical showers 



Type of 


Energy of 


Threshold 


Fraction of 


Number of 


Quality 


primary 


primary 


radius of 


showers 


showers 


factor 




(GeV) 


curvature 


accepted 


simulated 








(km) 


(%) 


(%) 




7— rays 


100 


11.7 


99 


200 


1.13 ± 0.13 


and protons 


250 




77.5 


200 




7— rays 


500 


8.8 


98 


200 


1.13 ± 0.13 


and protons 


1000 




75 


200 




7— rays 


1000 


8.2 


99 


100 


1.11 ± 0.18 


and protons 


2000 




80 


100 




7— rays 


1000 


6.3 


48 


100 


4.8 ± 2.6 


and Fe nuclei 


10000 




1 


100 




7— rays 


spectrum 


8.3 


100 


100 


1.13 ± 0.18 


and protons 


spectrum 




79 


100 





simulated in each case are also tabulated. The last row in the table shows the 
quality factor estimated for primaries chosen from a powerlaw spectrum as 
mentioned before. This is consistent with that for mono-energetic primaries 
showing that a change in the height of shower maximum with energy does not 
change the quality factor significantly. Fig. 1 shows the distributions of the 
fitted radii of curvature of the shower front for proton primaries and Fe pri- 
maries with respect to that for 7— ray primaries of equivalent Cerenkov yield. 
The vertical lines represent the threshold value with respect to which quality 
factors are estimated. One can see a rather high degree of overlap between 
the pair of distributions and hence the quality factors from this parameter are 
rather modest in value but are almost independent of primary energy. How- 
ever we will see in §7.1 that the use of a circular mask to limit the telescope 
opening angle improves the quality factor at all energies. Also the quality fac- 
tor against Fe primaries improves dramatically suggesting that this could be 
the ideal parameter to discriminate heavy primaries. 



5 Pulse shape parameters 
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Fig. 1. Distribution of the shower front radii for 7— rays (continuous line) and 
protons or Fe nuclei (dashed line) of respective energies as indicated in each panel. 
The last panel shows similar distributions for 7-rays & protons of primary energies 
chosen from a power law spectrum of slope -2.65 (see text for details). The vertical 
lines indicate the threshold values of parameters. 



5. 1 Observation level of Pachmarhi 



Even though pulse shape parameters contain information on the history of 
shower development in the atmosphere, their use for identifying the primary 
species was always in doubt especially at lower primary energies [27]. On 
the other hand, it was known that the muons generated in hadron initiated 
showers will reach the observation level several nanoseconds before the light 
from the electrons. Consequently, the light from muons would give rise to an 
unmistakable precursor which could be used as a discriminating parameter 
[29]. 

Here we tried to investigate the possibility of using pulse shape parameters 
to discriminate between 7— ray and hadron showers, using the same data 
set as was used for shower front paramemters. Distributions generated using 
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pulse shape parameters were derived from predicted lognormal distribution 
for individual detectors (see [11] for details). It may be recalled here that the 
decay time of the pulse need not be measured directly for this purpose. If 
one measures the mean arrival time and the RMS fluctuations in Cerenkov 
photon arrival times at a detector, it will enable us to derive the pulse shape 
parameters from the data directly [11]. Among the three pulse shape param- 
eters viz. the rise time, the decay time & the pulse width, decay time seems 
to be the most sensitive parameter which exhibits species sensitive behaviour 
at TeV energies. However, Cerenkov photon arrival time distribution yields 
one set of pulse shape parameters for each of the 357 detectors. So also the 
arrival time jitter (discussed in §6). Hence for the purpose of estimating the 
statistical errors on the quality factors, the sample size considered is the prod- 
uct of the number of detectors and the number of simulated showers. In order 
to take into account the core distance dependence of these parameters (see 
[11] for details) they are averaged over 160 consecutive detectors arranged in 
the order of increasing core distance. The resulting 22 (352/16, excluding the 
detectors at the shower core and the farthest 4) sets of parameters are treated 
as independent samples of the parameter under study. The distributions of 
the these samples of decay time both for 7— ray and proton primaries of var- 
ious energies are shown in fig. 2. The number of showers simulated are same 
as those listed in table 1. The distributions from hadronic primaries are in 
general characterized by longer tails compared to that from 7— ray primaries. 
There is a fair amount of overlap between the distributions of two types of 
primaries. All the relevant parameters including the fractions of 7— rays and 
protons whose decay times are less than the threshold value are listed in table 
2. In spite of the modest quality factor, one is able to reject nearly 96% of 
protons from the data while loosing nearly two thirds of the 7— ray signal. 

The other pulse shape parameters like the rise time and the pulse width (de- 
fined as the full width at half maximum, FWHM) for vertical showers have 
been found to be quite insensitive to the primary species and hence are not 
useful parameters to distinguish hadronic events. However as discussed in §7.3 
it can be seen that the pulse width could be a good parameter for inclined 
showers. 



6 Timing Jitter 

The Cerenkov light emitted at the top of the atmosphere reaches the observer 
later than that from deeper down. The latter appears at larger angles but 

2 The choice of this number is purely arbitrary. We chose this number because the 
spatial separation of the resulting 22 samples turned out to be reasonably uniform. 
However this grouping was not done or inclined showers. 
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Table 2 

Quality of pulse decay time as a discriminating parameter for vertical showers. The 
number of showers simulated are same as those indicated in table 1. 



Type of 


Energy of 


Threshold 


Fraction of 


Quality 


primary 


primary 


value 


showers 


factor 




(GeV) 


(ns) 


accepted (%) 




7— rays 


250 


3.2 


37.7 


1.67 ±0.02 


and protons 


500 




5.1 




7— rays 


500 


3.7 


35 


2.22 ± 0.03 


and protons 


1000 




2.5 




7— rays 


1000 


3.9 


30.4 


2.41 ± 0.06 


and protons 


2000 




1.6 




7— rays 


1000 


4.0 


32.1 


2.48 ± 0.06 


and Fe nuclei 


10000 




1.7 




7— rays 


spectrum 


3.9 


32.5 


1.57 ± 0.03 


and protons 


spectrum 




4.3 





arrives earlier because the particles travel faster than light in the atmosphere. 
It was suggested long ago that such small time differences could be useful 
in discriminating against proton showers [27]. We have seen before, that a 
majority of the Cerenkov photons originate at around the shower maximum 
as a result of which at observation level we observe the spherical front centred 
around that point. However the spread in the arrival time at any core distance 
is largely decided by the photons emitted elsewhere in the atmosphere. For 
example, at large core distances a bulk of the photons are emitted at lower 
atmospheric heights, largely from lower energy electrons undergoing multiple 
Coulomb scattering. In addition, photons seem to arrive at increasingly larger 
angles away from the core thus implying larger arrival times. In other words, 
the spread in the arrival times within a given shower at any core distance has 
the definite signature of the kinematics of the shower development. Hence we 
tried to search for species dependent signature in the arrival time spread of 
Cerenkov photons. We quantify the arrival time spread in terms of the RMS 
of the photon arrival times at a given detector. The ratio of this RMS to 
the mean arrival time, called the relative jitter, is used as a discriminating 
parameter to identify the primary species. 
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Fig. 2. Distribution of decay times for 7— rays (solid line) and protons or iron nuclei 
(dotted line). The vertical lines indicate the threshold values. 
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6.1 Core distance dependence 



It may be recalled that each telescope in the array considered here consists 
of compactly mounted 7 paraxial parabolic mirrors of diameter 0.9 m. The 
average arrival time of Cerenkov photons is calculated for each of the seven 
mirrors of a telescope, for a given shower. Then the mean arrival time for each 
telescope is calculated which is the average of seven individual averages. Time 
jitter is the RMS of seven averages for each telescope as mentioned before. 
Figure 3 shows a radial plot of the mean arrival times (ns, marked by +), time 
jitter (ns, marked by open triangles) and the relative jitter (which is the ratio 
of time jitter to the mean, asterisks) both for 1 TeV 7— rays (left) and 2 TeV 
protons (right). While there is a general increase in the jitter with increasing 
core distance for both the primaries, it shows a minimum at around the hump 
region for 7— ray primaries. This is expected since most of the photons here 
are emitted by the high energy electrons and also from a limited range of 
atmospheric heights. Such a signature is swamped by kinematic fluctuations 
in the case of hadronic primaries, in addition to exhibiting an increased time 
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Table 3 

Quality of jitter as a discriminating parameter. The number of showers simulated 
is 100 in all cases for each type of primary. 



Type of 


Energy of 


Threshold 


Fraction of 


Quality 


primary 


primary 


value 


showers 


factor 




(GeV) 




accepted (%) 




7— rays 


100 


0.23 


49.8 


2.83 ± 0.05 


and protons 


250 




3.1 




7— rays 


1000 


0.07 


67.5 


2.42 ± 0.03 


and protons 


2000 




7.8 




7— rays 


1000 


0.09 


91.3 


9.13± 0.25 


and Fe nuclei 


10000 




1.0 




7— rays 


spectrum 


0.08 


80.5 


1.85± 0.02 


and protons 


spectrum 




19 





spread [10]. 



6.2 Separation 0/7 - rays from proton primaries 

It may be seen from figure 3 that the relative time jitter for both 7— ray as 
well as proton primaries is almost independent of the core distance. Hence 
it is chosen as a parameter for 7— hadron separation even though it is not a 
necessary condition. 

Fig. 4 shows the comparative distributions of this parameter for 7— rays and 
protons (a & b) as well as 7— rays and Fe primaries (c) of equivalent Cerenkov 
yields. Also shown in the figure (d) are comparative distributions for 7— ray 
and proton primaries of varying energies selected from a power law distribution 
of slope -2.65. The energy bandwidths chosen for 7— rays and protons are 0.5 
-10 TeV and 1-20 TeV respectively. The results are based on 100 showers 
in each case for each species. 

It is well known that after shower maximum, 7— ray showers attenuate pro- 
gressively faster with atmospheric depth than do hadronic showers [36]. In 
addition, the lateral spread also is larger for hadron initiated showers. As a re- 
sult, the RMS spread in the time of arrival of photons at the observation level 
is expected to be larger for hadronic showers. Consequently, the distributions 
in fig 4 are pretty wide with a long tail for hadrons at all primary energies 
which makes it possible to discriminate them. 
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Fig. 3. Radial variation of the mean arrival time, ns (+), arrival time jitter, ns 
(open triangle) and the relative jitter (asterisk) for 1 TeV 7— ray (left) and 2 TeV 
proton (right) incident vertically at the top of the atmosphere. 

Table 3 shows the quality factors derived from relative jitter for vertical show- 
ers of 7— rays, protons and iron primaries. From the table it can be seen that 
using relative jitter as a discriminating parameter one can reject more than 
~80% of protons while retaining more than ~80% of 7— rays (from spectrum). 



7 Sensitivity to various operational parameters 



7. 1 Effect of telescope opening angle on quality factors 



The opening angle of a Cerenkov telescope is often limited by placing a circular 
mask at the focal point in front of the photocathode. This limits the arrival 
angle of the photons reaching the photocathode. In the absence of a mask the 
opening angle is limited by the photocathode diameter. In other words, the 
limiting mask is expected to reduce the mean arrival angle of photons, reduce 
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Fig. 4. Distributions of mean relative jitter for 7— rays (solid line) and protons 
or iron nuclei (dotted line) of different primary energies are shown, (d) shows the 
distribution for primaries with a power law (differential slope -2.65) distribution of 
energies in the range 500 GeV - 10 TeV for 7— rays and 1 TeV - 20 TeV for protons. 
The vertical lines represent the threshold values. 

the mean arrival time as well as increase the average production height. Table 4 
summarizes the effect on aforesaid shower parameters by placing a 3° (FWHM) 
mask at the focal point of a detector. This is to demonstrate qualitatively its 
effect on some of the parameters of shower development. Vertically incident 
showers are considered here. The changes in the mean angle of arrival and time 
is around 30% for 7— rays while it is significantly larger for proton primaries. 
The mean production height is much less sensitive to the presence of mask 
which increases the mean production height by about 10% and 14% for 7— ray 
and proton primaries respectively. 

The question we are addressing here is the effect of telescope opening angle 
on each of the quality factors discussed here. Tables 5, 6, 7 & 8 summarize the 
results relating to the role of mask when radius of curvature, pulse decay time, 
pulse width and relative jitter are used as parameters to discriminate against 
cosmic ray hadrons. The same data set is used to compare the results with 
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Table 4 

Effect of adding a mask at the focal point on some of the parameters for vertical 
showers 



Type of 


Energy of 


Mask 


Mean 


Mean 


Mean 


primary 


primary 


used 


arrival angle 


arrival time 


production 




(GeV) 


(FWHM) 


(degrees) 


(ns) 


height (km) 


7-ray 


100 


No mask 


1.3 ±0.27 


7.23 ±1.91 


11.04 ± 1.15 


7-ray 


100 


3.0° 


0.9 ±0.11 


5.11 ±0.66 


12.09 ± 1.13 


Proton 


250 


No mask 


1.64 ±0.42 


12.01 ± 18.29 


7.83 ± 1.11 


Proton 


250 


3.0° 


0.93 ±0.11 


7.11 ± 14.41 


8.91 ± 1.11 



and without mask. In the no-mask case all the Cerenkov photons incident on a 
detector, irrespective of their angle of incidence, are accepted. Pulse width was 
included in this case to show that if a mask is in use it could be a useful pa- 
rameter to discriminate against heavy primaries. Its sensitivity improves with 
increasing primary energy consistent with other such studies. From tables it 
can be readily seen that the use of a mask (5° FWHM) improves the efficiency 
of the radius of curvature as a discriminating parameter while the pulse shape 
parameters show a marked improvement only for heavy primaries. For the rel- 
ative timing jitter the efficiency as a discriminating parameter reduces when 
a mask is used. This, even though a marginal effect, is quite understandable. 
Limiting the photon arrival angles to near vertical direction reduces the RMS 
fluctuations because of limited range of pathlength differences. 



7.2 Altitude dependence of quality factors 

As the altitude of observation increases, the shower maximum for a given 
primary energy comes closer to the observation level. As a result, the lateral 
distribution of Cerenkov photons has been known to change with the altitude 
of the observation level. The core distance at which the hump appears as well 
as the prominence of the hump will be smaller with increasing altitude of 
observation [37,38]. Since the Cerenkov front is being intercepted at differ- 
ent observation levels during its propagation in the atmosphere, the shower 
parameters like the average arrival angle, time etc also will be different at dif- 
ferent observation levels. As we will see in §7.5 that the species discriminating 
efficiency of the parameter does vary, even though weakly, with core distance. 
Hence it would be interesting to see if their role depends on the observation 
level too. Therefore we studied the role played by the observation altitude in 
using the various types of parameters studied here. 
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Table 5 

Quality of radius of spherical wavefront as a discriminating parameter for vertical 
showers when 5° mask is used. 



Type of 


Energy of 


Threshold 


Fraction of 


Quality 


primary 


primary 


radius of 


showers 


factor 




(GeV) 


curvature 
(km) 


accepted 

(%) 




7— rays 


100 


12.5 


97 


1.37 ± 0.16 


and protons 


250 




50.5 




7— rays 


500 


9.7 


96.5 


1.35± 0.16 


and protons 


1000 




51 




7— rays 


1000 


8.8 


93 


1.33 ± 0.22 


and protons 


2000 




49 




7— rays 


1000 


8.7 


91 


9.1 ± 4.8 


and Fe nuclei 


10000 




1.0 




7— rays 


spectrum 


9.5 


100 


l.O ± u.z 


and protons 


spectrum 




61 




Table 6 

Quality of decay time as a discriminating parameter w 


ith 5° mask 


Type of 


Energy of 


Threshold 


Fraction of 


Quality 


primary 


primary 


value 


showers 


factor 




(GeV) 


(ns) 


accepted (% 


) 


7— rays 


250 


3.4 


52.8 


1.14± 0.01 


and protons 


500 




21.5 




7— rays 


500 


3.5 


44.9 


1.84 ± 0.02 


and protons 


1000 




6 




7— rays 


1000 


3.5 


40.1 


1.81±0.03 


and protons 


2000 




4.9 




7— rays 


1000 


4.8 


71.4 


6.98 ± 0.19 


and Fe nuclei 


10000 




1.1 




7— rays 


spectrum 


3.3 


35.2 


1.70 ±0.03 


and protons 


spectrum 




4.3 
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Table 7 

Quality of FWHM as a discriminating parameter when 5° mask is used. 



Type of 


Energy of 


Threshold 


Fraction of 


Quality 


primary 


primary 


value 


showers 


factor 




(GeV) 


(ns) 


accepted (%) 




7— rays 


250 


3.8 


98.2 


1.02±0.01 


and protons 


500 




92.5 




7— rays 


500 


2.0 


45.7 


1.08±0.01 


and protons 


1000 




18.1 




7— rays 


1000 


2.1 


39.3 


1.32±0.02 


and protons 


2000 




8.9 




7— rays 


1000 


2.0 


33.8 


3.31 ±0.09 


and Fe nuclei 


10000 




1.1 




7— rays 


spectrum 


2.2 


45 


1.07 ±0.01 


and protons 


spectrum 




17.6 




Table 8 










Quality of jitter 


as a discriminating parameter when a 5° 


mask is use< 


Type of 


Energy of 


Threshold 


Fraction of 


Quality 


primary 


primary 


value 


showers 


factor 




(GeV) 




accepted (%) 




7— rays 


100 


0.2 


58.2 


1.7 ± 0.02 


and protons 


250 




12.3 




7— rays 


1000 


0.04 


62.1 


2.1 ± 0.03 


and protons 


2000 




9.1 




7— rays 


1000 


0.05 


89.6 


8.8 ± 0.2 


and Fe nuclei 


10000 




1.1 




7— rays 


spectrum 


0.05 


86.4 


1.7 ± 0.02 


and protons 


spectrum 




26.7 





Tables 9, 10 & 11 summarize the quality factors for the pulse decay times as 
well as the relative timing jitter at two different observation levels viz. sea level 
and 2.2 km above mean sea level (for 500 GeV 7— rays and 1 TeV protons). 
In comparison with the quality factors for the same parameters in table 2 & 3 
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Table 9 

Quality of decay time as a parameter measured at sea level when used to discrimi- 
nate 500 GeV 7— rays from 1 TeV protons. 



Type of 


Mask 


Threshold 


Fraction of 


Quality 


primary 


used (°) 


value 


showers 


factor 




(FWHM) 


(ns) 


accepted (%) 




7— rays 


No mask 


3.2 


32.9 


2.82 ± 0.07 


and protons 


used 




1.4 




7— rays 


5° mask 


3.0 


37.5 


3.45±0.09 


and protons 


used 




1.2 




Table 10 

Quality of pulse decay time as measured at an altitude of 2.2 km 
discriminate 500 GeV 7— rays from 1 TeV protons. 


Type of 


Mask 


Threshold 


Fraction of 


Quality 


primary 


used (°) 


value 


showers 


factor 




(FWHM) 


(ns) 


accepted (%) 




7— rays 


No mask 


4.2 


31.5 


1.07±0.02 


and protons 


used 




8.6 




7— rays 


5° mask 


6.8 


97.3 


1.06 ±0.01 


and protons 


used 




84.1 





it can be seen that the quality factors from pulse decay time improve steadily 
with decreasing altitude while those from timing jitter are almost independent 
of altitude. In order to understand this we computed the quality factors, for 
detectors only around the hump region after taking into account the varying 
hump distances from the core at three altitudes. It is found that the quality 
factors so calculated for pulse decay time are respectively 8.1, 6.6, and 3.6 for 
sea level, 1 km and 2.2 km altitudes. As one would see later (table 16) that 
the quality factors computed exclusively for pre-hump & post-hump regions 
are poorer than that for the hump region. Hence the improvement in quality 
factor at lower altitudes is mainly due to the increased prominence of hump 
at lower atmospheric depths. 



The relative timing jitter, on the other hand seems to be much less sensi- 
tive to core distance (see figure 3) and hence doesn't vary significantly with 
observation altitude. 
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Table 11 

Quality of jitter as a discriminating parameter at two observation levels viz sea level 
& 2.2 km above mean sea level. Monoenergetic 7— rays of 500 GeV and protons of 
1 TeV incident vertically at the top of the atmosphere are considered. 



( 1 hQPT"\ratinvi 

V ./ Uot. 1 VojLIVJII 


1V_LOjOJ\. 


i 111 COllUlU. 




H T*5i r'tioTi of 

J. IdiUllUll Ul 


On n 1 


1 p vpl f] Yi nvp 

CliUU V c 


used (°) 


VcllllC 






fart or 

1CXV_, IU1 


nisi (km} 


fFWHM) 






(7, p) 




0.0 


None 


0.08 


52.5 


2.7 


3.18 ± 0.06 


1.07 


None 


0.09 


45.8 


1.1 


4.48 ± 0.12 


2.2 


None 


0.09 


41.2 


1.7 


3.16 ± 0.07 


0.0 


5.0 


0.05 


41.0 


3.3 


2.27 ± 0.03 


1.07 


5.0 


0.05 


48.9 


3.4 


2.67 ± 0.05 


2.2 


5.0 


0.05 


57.3 


4.5 


2.7 ± 0.04 



7.3 Incident angle dependence of quality factors 



At larger zenith angles rising absorption and the increasing distance from the 
shower maximum raise the energy threshold of the primary that are detected 
by an atmospheric Cerenkov telescope. After the shower maximum 7— ray in- 
duced showers attenuate progressively faster with atmospheric depth than do 
the hadronic showers. As a result one would expect a zenith angle dependence 
on the sensitivity of the parameters studied here. Qualitatively speaking in- 
clined showers at a given altitude behave similar to vertical showers at a 
lower altitude. As a result, one would expect the quality factor to improve 
with zenith angle, as it would with a decrease in the altitude of observation. 
The species sensitive imaging parameters like the azwidth for example, have 
been shown to be much less sensitive for primaries incident at angles > 30° 
[19]. Hence we wanted to check the possible dependence of the parameters, 
discussed here, on zenith angle of the primary. We have generated showers 
initiated by 7— rays and hadrons incident at the top of atmosphere at a fixed 
zenith angle of 30° to the vertical and uniformly distrbuted in azimuth. 

The quality factors have been derived at one test energy for comparison for 
all the parameters under consideration in the present study. The results are 
summarized in the following tables. 

Table 12 summarizes the efficiency of the radius of curvature of the light front 
as a discriminant for inclined showers of energy 500 GeV and 1 TeV for 7— ray 
and proton primaries respectively. The results are based on 200 showers each 
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Table 12 

Quality of radius of spherical wavefront as a discriminating parameter for inclined 
showers (incident at 30° to the vertical). 



Type of 


Energy of Threshold Fraction of 


Fraction of 


Quality 


primary 


primary radius of 7— rays 


protons 


factor 




(GeV) 


curvature accepted 
(km) (%) 


accepted 
(%) 




7— rays 


500 


10.9 94 


51 


1.32 ± 0.16 


and protons 


1000 








(without mask) 








7— rays 


500 


11.6 93.5 


46.5 


1.37 ± 0.16 


and protons 


1000 








(with 5° mas 


;k) 








Table 13 

Quality from decay time for showers incident at 30° to the vertical when used to 
discriminate 500 GeV 7— rays from 1 TeV protons. 


Type of 


Mask dia. 


Threshold Fraction of 


Fraction of 


Quality 


primary 


used 


value 7— rays 


protons 


factor 




(FWHM) 


(ns) accepted (%) 


accepted (%) 




7— rays 


No mask 


5.4 59.9 


11.3 


1.79 ± 0.02 


and protons 


used 








7— rays 


5.0° 


4.6 52.5 


8.2 


1.83 ± 0.02 


and protons 










of 7— rays and protons. It can be readily seen from a comparison with table 
1 that there is a significant improvement in the quality factors for inclined 
showers. This could be partly due to an increase in the range of fitted radii of 



proton showers. 

Tables 13 & 14 summarize the results of using two of the pulse shape pa- 
rameters viz. pulse decay time and pulse width for species identification for 
showers incident at 30° to the vertical. The results are based on 100 show- 
ers each of 7— rays & protons. Keeping in mind that the decay times are not 
grouped for inclined showers (unlike vertical showers) quality factor for decay 
time improves for inclined showers. 

Table 14 on the other hand has an interesting result. As mentioned in §5.1, 
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Table 14 

Quality of pulse width (FWHM) as a parameter from showers incident at 30° to 
the vertical when used to discriminate 500 GeV 7— rays from 1 TeV protons. 



Type of 


Mask dia. 


Threshold 


Fraction of 


Fraction of 


Quality 


primary 


used 


value 


7— rays 


protons 


factor 




(FWHM) 


(ns) 


accepted (%) 


accepted (%) 




7— rays 


No mask 


3.9 


84.3 


48.8 


1.21± 0.01 


and protons 


used 










7— rays 


5.0° 


3.4 


75.8 


33.6 


1.31 ± 0.01 


and protons 













pulse rise time and pulse width do not exhibit significant sensitivity to primary 
species for vertical showers. However for inclined showers pulse width becomes 
a useful parameter as shown in table 14. This is consistent with the results 
of Roberts et al. [28], who showed that the quality factor for pulse width 
peaks at around an inclination of ~ 35° to the vertical at a given altitude of 
observation. 

The relative shower to shower fluctuations in pulse shape parameters have 
been found to be less for inclined showers [11] unlike at higher energies where 
the trend reverses [28]. Table 15 summarizes the quality factors using the 
timing jitter as the parameter. It may be seen from a comparison of tables 3 
& 15 that the species sensitivity of the timing jitter falls with zenith angle of 
the primary. 

The sensitivity of each of the parameters for inclined showers has been further 
studied after introducing a focal point mask. The results are again summa- 
rized in tables 12-15. The opening angles of inclined showers are smaller than 
that for vertical showers because, inclined showers reach shower maximum at a 
higher altitude where the Cerenkov angle is smaller. Further, the lower energy 
electrons propagate in a lower density medium thus undergoing less frequent 
Coulomb scattering. Consequently, the Cerenkov emission from inclined show- 
ers is more collimated compared to that from a same energy primary in vertical 
direction. Hence the use of a focal point mask is not expected to change the 
quality factors significantly since a smaller fraction of photons are obstructed 
by it. The results in tables 12-15 show no major change in quality factors when 
focal point masks are introduced, essentially supporting the above argument. 

For reasons mentioned above, the differential pathlengths of Cerenkov photons 
also reduce resulting in a poor sensitivity for timing jitter for inclined showers 
compared to that for vertical showers. The results of table 15 support this 
conclusion. Once again the use of a focal point mask has no significant effect 
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Table 15 

Quality of relative jitter as a discriminating parameter for 7— ray (500 GeV) and 
proton (1 TeV) primaries incident at 30° with respect to the vertical. 



Type of 


Mask dia. 


Threshold 


Fraction of 


Fraction of 


Quality 


primary 


used 


value 


7— rays 


protons 


factor 




(if any) 




accepted (%) 


accepted (%) 




7— rays 


No mask 


0.006 


66.3 


37.5 


1.08±0.01 


and protons 


used 










7— rays 


5.0° 


0.009 


99.8 


98.7 


1.01±0.01 


and protons 













on the quality factor. However at higher primary energies it is observed that 
the spread in the arrival times of Cerenkov photons increases with the mass 
of the primary and the zenith angle [28]. This is due to the increased muon 
production in cascades generated by the heavier primaries. 



7.4 Species dependence of quality factors 

Each of the three species sensitive parameters under study here were applied to 
Fe primaries as well. It can be easily seen that in all the cases these parameters 
are more efficient in discriminating heavy primaries as compared to protons. 
The radius of the shower front, for example, shows a higher sensitivity for 
Fe primaries which is a reflection of the fact that, an Fe shower reaches its 
maximum development at a higher altitude compared to that of a lower Z 
primary. The efficiency also improves significantly with the use of a mask, a 
feature similar to that of proton primaries. 

The pulse shape parameters on the other hand do not exhibit a large change 
in sensitivity for high Z primaries. However the use of mask shows a drastic 
improvement in sensitivity of decay time as the quality factor changes from 
~ 2.5 (with out mask) to ~ 7 (with a 5° mask). 

The photon arrival time jitter is particularly sensitive to high Z primaries. 
This perhaps is quite understandable since the Fe shower if considered as 
superposition of proton showers is expected to be more spread out laterally 
and hence the timing jitter will deviate more than that for 7— ray showers. 

This demonstrates that the parameters discussed here are sensitive not only to 
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Table 16 

Quality of decay time as a discriminating parameter for 7— ray (500 GeV) and 
proton (1 TeV) primaries at three different core distance ranges: pre- hump, hump 
& post hump 





TVTa^lc Hi a 


r TVi rpc Vl ol H 


T^rartion of 


T^rartion of 


Onalitv 




used 


Veil 116 


/ icy o 


nrot.on s 


factor 

-LCI V_< L Ul 




(if anv) 


(us) 


accepted (%) 


accepted (%) 




Pre-hump 


No Mask 


4.3 


88.5 


30.5 


1.60 ± 0.05 


Hump 


No Mask 


3.3 


80.5 


1.5 


6.57 ± 0.51 


Post-hump 


No Mask 


8.4 


53.0 


14.0 


1.42 ± 0.06 


Pre-hump 


5.0° 


3.8 


73.0 


19.5 


1.65 ± 0.06 


Hump 


5.0° 


2.8 


91.5 


1.0 


9.15 ± 0.85 


Post-hump 


5.0° 


6.0 


93.5 


41.5 


1.45 ± 0.04 



proton primaries but also to high Z primaries, unlike Hillas's image parameters 
which are not very sensitive to heavy primaries of cosmic rays. 



7. 5 Core distance dependence of quality factors 

There are mainly three effects that determine the lateral distribution of Cerenkov 
light at the observation level: the finite Cerenkov angle, Coulomb scattering of 
the progenitor electrons and the transverse momentum p t , in hadronic inter- 
actions. While the first effect is mainly responsible for the proverbial 'hump' 
in the lateral distribution of 7— ray initiated showers the third effect is respon- 
sible for the observable differences between the two types of primaries and the 
second effect is responsible for smoothing out the other two effects. As a re- 
sult, one expects species specific parameters to show a dependence on the core 
distance of the detectors which measure these parameters. Table 16 summa- 
rizes the relative efficiencies of one of the parameters, viz pulse decay time as 
measured by detectors at pre-hump, hump & post-hump regions of showers. It 
can be seen that the difference between 7— ray and proton initiated showers is 
maximum at around the hump region resulting in maximum quality factor at 
a given primary energy. At the other two regions, viz. pre-hump & post-hump 
regions, the quality factors are comparable. However in practice it may not be 
always possible to detect the hump. Hence it may be advisable to limit the 
core distance of the triggering events to within the hump region at a given 
altitude of observation which is quite feasible since hump radius is not very 
sensitive to primary energy. 
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8 Discussions 



In the present study we did not take into account the effect of detector re- 
sponse on the parameters selected here. It is implicit that in an experiment 
the pulse shapes have to be measured using high bandwidth electronics to 
minimize the shape distortion due to instrumental response so that the pulse 
shape parameters could be effective discriminators. However while applying 
these techniques to real data, it is inevitable that the theoretical parameters 
discussed here be suitably corrected for instrumental response used in the sys- 
tem. The threshold values of the parameters have to be derived for a given 
detection system and then applied to the data in order to reject events of 
hadronic origin. 

The quality factors derived from the current studies have been tested for 
stability by estimating the same for subsets of the totality of data used. Hence 
their values are not critically dependent on the sample size. However it must 
be mentioned we chose to use that value of the parameter averaged over 16 
detectors, thus reducing the 357 radial samples to 22. Obviously, the quality 
factors will improve if one averages over a larger number of detectors, as has 
been seen in our analysis. This only means that arrays with larger number of 
detectors can reject hadronic showers with a better efficiency. 

The effect of mask has been studied qualitatively here. In order to apply the 
technique to real data one has to use the actual mask diameter. It is possi- 
ble, in principle, to optimize the mask diameter for the best hadron rejection 
efficiency 

8. 1 Radius of curvature of the shower front: 

From the table 1 it can be seen that the radius of curvature of the Cerenkov 
front can enable us to reject around 20% of proton events at all energies. The 
rejection efficiency could improve to 50% if one uses a mask to limit the photon 
incidence angle. On the other hand this parameter is very sensitive to heavy 
primaries and hence one can reject up to ~ 99% (with or without a 5° mask) 
of them. 

At higher primary energies (> 20 TeV), the height of shower maximum has 
been shown to be directly proportional to the slope of the Cerenkov photon 
lateral distribution [39]. This has been suggested to be a good parameter 
to distinguish between 7— rays and hadrons. The present results show that 
the height of shower maximum, measured as the fitted curvature of the light 
front, is a reasonably good parameter even at lower energies and the sensitivity 
improves for inclined showers. As already mentioned, it is particularly useful 
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to discriminate against heavy primaries. It can be readily seen from tables 1 
& 5 that the quality factor is sensitive to photons arriving at large incidence 
angles. This is particularly so for heavy primaries. 

However the conventional method of measuring the height of shower maximum 
at higher primary energies is by measuring the pulse width [40]. This technique 
seems to work only at large core distances and at primary energies > 10 15 eV 
as the two parameters do not show any correlation at primary energies studied 
here. 

8.2 Pulse shape parameters 

There are not many results published in the literature regarding the useful- 
ness of the pulse shape parameters to discriminate hadronic showers from 
electromagnetic showers at lower primary energies. Rodriguez-Fries et al. [41] 
reported that the pulse width could be used as a discriminating parameter at 
primary energies > 100 TeV and their results show that sensitivity of this pa- 
rameter is too poor at 10 TeV. This conclusion is consistent with the present 
results. 

Patterson & Hillas [27] too have concluded that the pulse widths of simulated 
proton showers are not very different from those for electromagnetic showers 
and hence unlikely to be useful as a discriminating parameter consistent with 
the conclusions of the present study. However we find that it could be a useful 
parameter for distinguishing inclined hadronic showers. 

The present results show that the pulse width and pulse decay time to be 
good discriminating parameters when used with a mask. However our study 
shows no sensitivity to pulse rise time in contrast to the conclusions of Roberts 
et al. [28]. Roberts et al. on the other hand do not use fall time to achieve 
discrimination against hadrons for they find that optimum quality factor was 
achieved only by rejecting an excessive fraction of 7— ray cascades, unlike in 
the present studies. They also find that rise time cut discriminates best against 
heavy primaries while we find that the decay time discriminates best against 
heavy primaries. The 7— ray acceptance fractions when decay time is used as 
a discriminating parameter is in the range 30-38% in contrast to the results 
of Roberts et al. [28]. 

8. 3 Timing jitter 

Cabot et al. [29] have estimated the relative mean times of arrival of Cerenkov 
photons at various core distances. The present results shown in fig. 3 (shown 
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Fig. 5. Variation of mean relative jitter for 7— rays (left) and protons (right) as 
a function of primary energy. The primary energies were selected randomly from 
power law distribution (differential slope -2.65). The energy band widths are 500 
GeV - 10 TeV for 7— rays and 500 GeV - 20 TeV for protons. The fitted slopes are 
indicated for each species. 



as +) agree well with theirs. The mean arrival times are fairly independent of 
the primary species. However the RMS fluctuations in arrival times depend 
on the range of differential pathlengths which in turn depend on the details 
of the interaction kinematics. Hence they bear the signature of the primary 
radiation. 

One can ask the question how does the relative jitter vary with primary energy. 
It could be seen from figures 4a, 4b & 4d that the distributions of relative jitter 
become narrower with increasing energy for both 7— ray and proton primaries 
while average value of the relative jitter falls with increasing primary energy 
far more rapidly for proton primaries than for 7— rays. Figure 5 shows the 
primary energy dependence of relative jitter averaged over all core distances. 
Hence the quality factor for 7— ray and proton primaries is expected to fall 
with increasing primary energy as could be seen from table 3. 
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In an experiment it is not readily obvious how to measure the mean arrival time 
and the timing jitter. Normally most of the experiments use discriminators 
which trigger at known thresholds. Knowing that a typical Cerenkov pulse 
can be well described by a lognormal function (modulated by the phototube 
response) it is possible to estimate its coefficients knowing when the pulses 
from individual elements cross the preset threshold as well as the pulse height. 
The coefficients of the lognormal functions are related to the mean arrival 
time and its RMS value. It is therefore essential to have independent multiple 
sampling of the threshold crossing time and photon density from a single 
Cerenkov telescope. Further details are beyond the scope of this paper. 



9 Conclusions 

Three different measurable parameters based on Cerenkov photon arrival times 
at detectors at different core distances at the observational level have been 
found to be useful in distinguishing between electromagnetic and hadronic 
showers. The sensitivity of these parameters peaks at around the hump region 
and hence best discrimination would be obtained by limiting the measurements 
from detectors around the hump region. The sensitivity of these parameters 
also seems to increase with the zenith angle of the primary at the top of the 
atmosphere. This is a useful property in contrast to the imaging technique 
which is most sensitive at small zenith angles [19]. Similarly, the use of a cir- 
cular mask at the focal point also increases the quality factor in some cases, 
especially for the curvature of the shower front. However the last two conclu- 
sions are not true to relative timing jitter in which case the sensitivity falls 
with increasing angle of incidence with respect to the vertical and marginally 
by the use of a mask. 

Although the quality factors for individual parameters are not very large by 
themselves, by applying them in tandem would greatly augment the detec- 
tion sensitivity of ground based VHE 7— ray telescopes designed to exploit 
the wavefront sampling techniques. Table 17 demonstrates the improvement 
in the quality factor when the decay time and the timing jitter are applied 
successively to the same sample of showers. The dramatic improvement in the 
quality factor exhibits the orthogonal nature of the parameters. 

The separation efficiency seems to decrease with increasing altitude of obser- 
vation mainly because of the decreasing prominence of the hump in the case 
of electromagnetic showers. This, once again supports the earlier conclusion 
that the Cerenkov photons around the hump region are more sensitive to pho- 
tonic primaries. Hence sampling of photons from this region of the light pool 
is recommended for better discrimination of hadronic primaries. 
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Table 17 

Table showing the improvement in quality factor when two parameters (viz. photon 
timing jitter and pulse decay time) are used in tandem to discriminate against 
hadronic events. For this 100 showers each of 7— rays and protons of energy 500 
GeV and 1 TeV respectively as sampled at sea level are used. 



Parameter 


Threshold 


Fraction of 


Fraction of 


Quality 


Type 


value 


7— rays 


protons 


factor 






accepted (%) 


accepted (%) 




Pulse Decay 


3.2 ns 


32.9 


1.4 


2.82±0.07 


Time 










Photon Timing 


0.08 


52.5 


2.7 


3.18±0.06 


Jitter 










Decay Time 


3.2 ns 


26.8 


0.05 


12.6±1.6 


& Timing Jitter 


0.08 
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